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r^- Bianchi VII/i models have been recently proposed to explain potential anomalies in the CMB 



< 



o 



anisotropy as observed by WMAP. We investigate the violation of statistical isotropy due to an 



^\i ! embedded Bianchi VII/i templates in the CMB anisotropy maps to determine whether the existence 

> 

Q^ ' of a hidden Bianchi template in the WMAP data is consistent with the previous null detection of 

!>■ 

CN ' the bipolar power spectrum in the WMAP first year maps. We argue that although correcting the 

^^ , WMAP maps for the Bianchi template may explain some features in the WMAP data it may cause 

>-^ ' other anomalies such as preferred directions leading to detectable levels of violation of statistical 

0^' isotropy in the Bianchi corrected maps. We compute the bipolar power spectrum for the low 
density Bianchi VII/j models embedded in the background CMB anisotropy maps with the power 



j^ [ spectrum that have been shown in recent literature to best fit the first year WMAP data. By 

^^ examining statistical isotropy of these maps, we put a limit of {jt)o < 2.77 x 10~"'^'^(99%CL) on 

j_j ■ the shear parameter in Bianchi Yllh models. 
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I. INTRODUCTION 

Friedmann- Robert son- Walker (FRW) models are the simplest models of the expanding 
Universe which are spatially homogeneous and isotropic. When proposed, the principal 
justification for studying these models was their mathematical simplicity and tractability 
rather than observational evidence. However, we now have observational evidence from 
the isotropy of CMB and large scale structure that the universe on large scales must be 
very close to that of a FRW model. However, there still remains some freedom to choose 
homogeneous models that are initially anisotropic but become more isotropic as the time goes 
on, and asymptotically tend to a FRW model. Bianchi models provide a generic description 
of homogeneous anisotropic cosmologies and they are classified into 10 equivalence classes 
[l6|. Among these models it is reasonable to consider only those types that encompass 
FRW models. These are types I and VHq in the case oi k = 0, V and VII^ in the case of 
k = —1, and IX in the case oi k = +1. Bianchi models which do not admit FRW solutions 
become highly anisotropic at late times. Type IX models re-collapse after a finite time and 
hence do not approach arbitrarily near to isotropy. Also models of type VII^ will not in 
general approach isotropy ll(]. The most general Bianchi types that admit FRW models are 
Bianchi types Vllh and IX. These two types contain types I, V, VHq as special sub-cases. 
An interesting feature of these models is that they resemble a universe with a vorticity. 
It is interesting to determine bounds on universal rotation from cosmological observations 
because the absence of such a rotation is a prediction of most models of the early universe, 
in particular, within the paradigm of infiation. 

CMB anisotropy is a powerful tool to study the evolution of vorticity in the universe 
because models with vorticity have clear signatures in the CMB. The vortex patterns that 
are imprinted by the unperturbed anisotropic expansion are roughly constructed out of two 
parts: production of pure quadrupole variations, or focusing of the quadrupole pattern in 
open models and a spiral pattern that are the characteristic signatures of VHq and Vllh 
models respectively l|, [uJ, [l5|. There are distinct features in each Bianchi pattern. For 



example the temperature anisotropy pattern for Bianchi Vllh universes is of the form 

—^(^,0) = /i(^)sin0 + /2(^)cos0. (1) 

J-o 

The functions fi{0) and /2(6') depend on the parameters of the particular Bianchi model and 
must be computed numerically Ij. In a pioneering work, analytical arguments were used 



to find upper bounds on the amount of shear and vorticity in the universe today, from the 
absence of any detected CMB anisotropy [ll|. A detailed numerical analysis of such models 
used experimental limits on the dipole and quadrupole to refine limits on universal rotation 
PJ. After the first detection of CMB anisotropy by COBE-DMR, Bianchi models were again 
studied by fitting the full spiral pattern from models with global rotation to the 4-year DMR 
data to constrain the allowed parameters of a Bianchi model of type Vllh [Sl, l39(]. Recently 
Bianchi Vllh models were compared to the first-year WMAP data on large scales and it was 
shown that the best fit Bianchi model corresponds to a highly hyperbolic model {Qq = 0.5) 
with a right-handed vorticity (■^)o = 4.3 x 10^^° 35|. They also found that 'correcting' 
the first-year WMAP data including the ILC map|61l] for the Bianchi template, makes 
the reported anomalies in the WMAP data such as alignment of quadrupole and octupole 
disappear and also ameliorates the problem of the low observed value of quadrupole [62 1. 
In a more recent work using template fitting it was found that although the "template" 
detections are not statistically significant they do correct the above anomalies 4J]. On the 
other hand, it was also found by [35] that deviations from Gaussianity in the kurtosis of 
spherical Mexican hat wavelet coefficients of the WMAP first year data are eliminated once 
the data is corrected for the Bianchi template. In a recent work the effect of this Bianchi 
correction on the detections of non- Gaussianity in the WMAP data was investigated [4g] and 
was shown that previous detections of non- Gaussianity observed in the skewness of spherical 
wavelet coefficients reported in [43] are not reduced by the Bianchi correction and remain at 
a significant level. It has been argued in Ref. [37| that increasing the scaling of the template 
by a factor of 1.2 makes the spot vanish. Although a stronger pattern may eliminate the 
non-Gaussian spot, it would tend to make the resultant map more anisotropic. Therefore 
there is a limit to the strength of anisotropic patterns that can be hidden in the WMAP 
data. 
teams 

In this paper we show that the statistical isotropy violation of a hidden pattern in the 
CMB map, such as that of Bianchi universe, can be quantified in terms of the Bipolar power 
spectrum. More specifically, we test the consistency of existence of a hidden Bianchi template 
in the WMAP data proposed in recent literature against our null detection of bipolar power 



?he above analysis was redone after the release of 3 year WMAP data 
71, l3^ and the previous conclusions were confirmed. 



33], by two 



spectrum (BiPS) in the WMAP first year and three year maps 26 



27 



281 ]. The bipolar 



power spectrum is a measure of statistical isotropy in CMB anisotropy maps and is zero 



when statistical isotropy obtains. Properties of BiPS have been studied in great details in 
[30| and [2|. We show that although correcting the WMAP first year maps for the Bianchi 
template may explain some features in the WMAP data63||, this is done at the expense 



of introducing some anomalies such as preferred directions and the violation of statistical 
isotropy into the Bianchi corrected maps. This violation is stronger in the case of enhanced 
Bianchi templates proposed by [STJ. By testing statistical isotropy of Bianchi embedded 



CMB maps, we put a limit of (f )o < 2.77 x 10-^^{99%CL) on Bianchi Vlh models. ¥. 
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a is the shear and H is the Hubble constant. Our result is marginally consistent with 

n 

but not with the enhanced Bianchi templates of [37|. 

The rest of this paper is organized as the following. Section 2 is a brief review of Bianchi 
classification of homogeneous spaces. Section 3 is dedicated to deriving the CMB patterns 
of Bianchi Vllh models. In § 4 statistical isotropy of Bianchi Vllh models are studied and 
compared to the null bipolar power spectrum of the WMAP data. Finally in § 5 we draw 
our conclusion. 

II. BIANCHI MODELS 

Today, the Bianchi classification of homogeneous spaces is based on a simple scheme for 
classifying the equivalence classes of 3-dimensional Lie algebras [6^]. This scheme uses the 
irreducible parts of the structure constant tensor under linear transformations. Following 
[l6|, we decompose the spatial commutation structure constants C^^ into a tensor, n"*, and 
a vector, ab 

CI, = edbcn"" + (5>, - 5ta, (2) 

where eabc is the 3-dimensional antisymmetric tensor and n""^ and a^ are defined as 

«6 = l^Cba (3) 

ab _ l^(a Med 

The structure constants C^, expressed in this way, satisfy the first Jacobi identity. Second 
Jacobi identity, shows that a^ must have zero contraction with the symmetric 2-tensor n"^; 

Ct^bCU = =^ ^^'^^ = 0- (4) 



We choose a convenient basis to diagonalize n'^* to attain n°-^ = diag(ni, n2,n3) and to set 
af, = (a, 0, 0). The Jacobi identities are then simply equivalent to riia = 0. Consequently we 
can define two major classes of structure constants 
Class A : a = 0, 
Class B : a 7^ 0. 

One can then classify further by the sign of the eigenvalues of n*^* (signs of ni, n2 and 
n^). Parameter h in class B, where a is non-zero, is defined by the scalar constant of 
proportionality in the following relation 

abac = -eukecjin'^ n^^ ■ (5) 

In the case of diagonal n""^, the h factor has a simple form h = a^ /{n2n^)- 

III. CMB PATTERNS IN TYPE VII;, BIANCHI MODELS 

This is the most general family of models which includes the k = —1 FRW solutions. It 
has some of the features of both types V and type VIIq- There is an adjustable parameter h 



in this model which is given by the square o 
model, the expansion scale factor, a, and 13 



" structure constants. In an unperturbed FRW 
65|] are given by [47[, 



e 



a 



/ll/2fi0 . ^2 fh^'^r 



sinhM^— , /5 = 0. (6) 



We introduce a factor x = HQe°'° which is related to the parameter h by 



As we said before, x has no physical effect on the FRW models. It can be seen from the 
fact that the present value of the scale factor can be arbitrarily chosen, but from eqn. ([6]) 
we have 

e"« = ^, (8) 

and hence, we see that parameter x can be scaled out of the solution. However, for large 
values of x (or equally h) , the models are similar to those of type V. As the present density 
tends towards the critical density, fio ^ 1, and h —>■ in such a way that x remains finite, 
the behavior of the models tend to that of type VIIq. Behavior of these models with different 
parameters can be seen in figure [H 



Null geodesies have a complicated form and are given by 



tan 



tan I -^ I exp [-(r - ro)v^] 

+ {t-tq) -= In {sin^ f y I + cos" 



(9) 



exp[2(r-ro)v^]}. 



They start near the South Pole and spiral in the negative 0-direction up towards the equator, 
and then spiral in the positive 0-direction up towards the North Pole. Small anisotropies 
can be treated as perturbations to FRW model. The value of d has been calculated for 

n ri 

them and is given in [Hj. Here we follow |l| who used these to calculate the contribution 
to temperature anisotropies from the vorticity alone. They argue that the inclusion of 
other pure shear distortions which are independent of the rotation could only make the 
temperature anisotropies larger. And hence their results would give the maximum level of 
vorticity consistent with a given value of temperature anisotropy. 

In type VII/i, there are two independent vorticity components uj2 and uo^. They are given 
in terms of the off-diagonal shear elements as 

(3/l - l)ori3 - 4/ll/Vi2 



OJ2 



^3 



:i - 3/i)ai2 - 4/iVVi3 



(10) 



The observables in this model are two dimensionless amplitudes: (^)q and (^)q- The 
vorticity is given by 

(11) 



uo 



]-e'%l + hf'^[{^ 



\2ll/2 



IM2J- + M']" 

where U2 and M3 are the velocity components and their present values are given by 



("^)° 3x^0 


r/'^i2\ , 0,1/2 /^i3\ ■ 
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the present value of vorticity is 


\HJo 6x2 


1 + 9/1)1/2 


Y^12\2 /( 

\H Jo V 


Tl3\2" 

H )o_ 



(12) 



1/2 



(13) 



To first order, temperature fluctuations are given by 

AT(6'o,0o) / irr\ / iTT\ 



pVVifcdt, 



(14) 



where 6*0 and 0o are related to the actual observing angles by 

9 = n - 9o, (f) = TT + (f)Q 



(15) 



Substituting for null geodesies from eqn. (^ and for velocities from eqn. (IT^ into eqn. (IT^ . 
we will obtain 

AT{9oAo) 



+ 






cos 00, 



(16) 



where coefficients ^(6*0) and -8(6*0) are defined by 



^(^0) = Ci [sin 00 -C2 (cos t/^s- 3/1^/2 sin Vs)] 
s(l — s^) sin0dr 



(17) 



+ Q 



Bi 



^oj 



l + s2)2sinh^/iV2r/2' 
Ci[3/i^/2 sin 00 - C2(sinV^s + S/i^/^ cos^^^)] 
p s(l-s2)cos0dr 



,,^ (l + s2)2sinh^/iV2r/2' 
The limits of integration are defined as 



ro 



= 2/^-^/2 gjni^-i (^-1 _ 1) 



1/2 



(18) 



te 



2/i"^/2 sinh 



-1 / ^^0 



and constants Ci, C2, C3, s and ■?/' are defined by 



1 + Ze 



Ci = (3fioa;)-^ 
02 — 



(19) 



1 + 4 



'9' 



s = tan I - 



(y lexp[-(r-ro)v^], 
^ = 0o + (r-ro)--^ln{sin2('|')+cos2('|')exp[2(r-ro)v^]}. 
The expression for — LgiW ^^^^ ^^ written in a compact form 



^{9o,<Po) = iA' + By/'(^) cos (00 + 
Jq Vii/o 



where is defined as 



cos I 



{^) B - {^) a] (A' + B')-'/' 



(20) 



(21) 
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FIG. 1: Bianchi patterns in type VII/i models for x = 0.1, 0.5, 1, 5 (from top to bottom) and 
r^o = 0.1, 0.5 and 0.9 (from left to right). In low density models, the pattern is mostly focused in 
one hemisphere but for nearly flat models the spiral patterns widen and extend to both hemispheres. 



and a is given by 



a = a 



12 



o 



13- 



(22) 



Eqn. (1201) helps to understand the behavior of the CMB pattern in these models. If we look 
around any circle at a given 9q on the sky, the temperature variation will have a pure cos 0o 
behavior, cj) determines the relative orientation of adjacent ^o =constant rings. On the other 
hand, B{9q) determines the amplitude of ^ and the focusing into a hot spot [l[. Some of 
these maps have been shown in figure [H 



IV. UNVEILING HIDDEN PATTERNS OF BIANCHI VII 



h 



Recently Bianchi VII/i models were compared to the first-year WMAP data on large scales 
and it was shown that the best fit Bianchi model corresponds to a highly hyperbolic model 
{VLq = 0.5) with a right-handed vorticity (^)o = 4.3 x 10^^° [35]. There are several aspects 
to this hidden pattern which should be mentioned. First of all, the proposed model is an 



8 



extremely hyperbolic model and does not agree with the location of the first peak in the 
best fit power spectrum of the CMB. In the regime of Vllh models, one should study the 
nearly flat models in order to have a consistent angular power spectrum. But the problem 
is that as it has been shown in Fig. [T], the spiral patterns in high density models are not 
concentrated in one hemisphere and can't be used to cure the large-scale power asymmetry 
in the WMAP data reported by [1^|. 

Second is that the Bianchi patterns we studied so far, are only valid for a matter domi- 
nated universe. Patterns must be recalculated in a universe with a dark energy component if 



one wants to compare 
dominated universe 



;hem with the observed CMB data which is believed to be there in a A 
. These issues are addressed in great details by [36] and they conclude 
that the "best-fit Bianchi type Vllh model is not compatible with measured cosmological 
parameters" . 

For these reasons, we study the type Vllh Bianchi models only as hidden patterns in the 



CMB anisotropy maps. We pay our attention to the specific model proposed by [35|| and in 
the next two sections, we address the question whether and to what extent one would be 
able to discover this pattern and other hidden anisotropic patterns in the CMB anisotropy 
maps. We carry out our analysis on full sky CMB maps. It has been shown in our previous 
papers 



30], that for a masked CMB sky maps, the Bipolar power spectrum has a calculable 
specific form. In that case, the violation of SI is measured with respect to the 'bias' that 
arises for the masked sky map. The excess BiPS signature would certainly depend to some 
extent on the specific orientation of the Bianchi template in the sky, given by the extent 
to which the pattern is covered by the mask. In this work we choose to keep our analysis 
independent of the orientation of the hidden pattern. 



Bipolar Power Spectrum Analysis 

We choose a Bianchi Vllh model with (^)q = (^)o ~ (l')o' "^^^ temperature fluctu- 
ations induced by this model are given by eqn. (IT6l) and will be 

AT(^o,0o; 



{j^)jMOo) + fi(^o)]sin0o (23) 



To 

+ [B{eo) - A{eo)]cos<Po, 

i^ower spectrum of the above Bianchi-induced fluctuations can be computed analytically 
l|, |46| and is given by 

where 



/ 2/ + 1 r^ 

'' = ^i^i{rTT)L '^(^)^'(-^'')-'"'^ (25) 

/ 2/ + 1 r 

The majority of the power in this Bianchi template is contained in multipoles below / ~ 20. 
Therefor to study this particular model we do not need high resolution CMB anisotropy 
maps. 

To do a statistical study of the Bianchi patterns, we generate the pattern for a given 
model. This pattern is given by the shear components, x parameter and the VLq. We will 
then simulate random CMB maps from the best fit Ci of the WMAP data 5J]. We add 
these two maps with a strength factor a which will let us control the relative strength of 
the pattern and the random map (see Fig. [2]). The resultant map is then given by 




+ M 



FIG. 2: Adding a pattern template with a strength a to a random realization of the CMB anisotropy 
map. This introduces a preferred direction in the map and violates the statistical isotropy. Maps 
are rotated to place the Galactic pole at the center for clear illustration. 

AT(n) = AT^*^^(n) + aAT^'''"'='^*(n), (26) 

10 



and hence the power spectrum of this map will be given by 



I = Ci + a Ci 



a is related to (f )o through (f )o =^ where Tq = 2.73 x lOV^ 



(27) 




'f.m. 
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FIG. 3: Four Bianchi added CMB maps with different strength factors, rotated to place the Galactic 
pole at the center for clear illustration. The Bianchi template in the above maps has been computed 
with X = 0.55, ilo = 0.5 and strength factors {§)o = 1.83 x 10"^ {top left), {jj)o = 1.09 x 10"^ [top 
right), {jj)o = 7.3 x 10"^° {bottom left), (-|)o = 3.66 x 10"^° {bottom right). Ahhough the pattern 
can hardly be seen in the fourth map, it has a considerable non-zero bipolar power spectrum (see 
Fig ED. 

We compute the bipolar power spectrum (BiPS) for the Bianchi added CMB anisotropy 
maps generated above. The unbiased estimator of BiPS is given by, 



Hi 



ll'M 



WiWv\Y,^i^^v^'Cw'? -^^. 



(28) 



where Wi is the window function in harmonic space, B^ is the bias that arises from the SI 
part of the map and is given by the angular power spectrum, Ci and Cf^^ are Clebsch-Gordan 
coefficients. Bipolar powe r sp ectrum is zero for statistically isotropic maps and has been 
studied in great details in 30| and |2]. Although BiPS is quartic in aim, it is designed to de- 



tect SI violation and not non-Gaussianity [25|, 



30 



31| . Since BiPS is orientation independent 
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[30|, we do not need to worry about the relative orientation between the background map 
and the Bianchi template. For different strength factors we simulate 100 Bianchi added 
CMB maps|67| at HEALPix resolution of Ngi^e = 32 corresponding to Imax = 95. Some 
of these maps are shown in Figure [31 BiPS of each map is obtained from the total aim 



which according to eqn. 



are ai^ 



-.CMB 
''Im 



I ^ „ Bianchi 



. We compute the BiPS for each 



map using the estimator given in eqn. (1281) . We average these 100 BiPS and compute the 
dispersion in them. The dispersion is an estimate of la error bars. At last we use the 
total angular power spectrum given in eqn. (pTj) to estimate the bias. We use the best 
fit theoretical power spectrum from the WMAP analysis, M], as Cf^^^ in computing the 
bias. Some results o " 
been discussed in |25|, 

order to concentrate on a particular /-range. This is proved to be a useful and strong tool 
to detect deviations from SI using BiPS method. In other words, multipole space windows 
that weigh down the contribution from the SI region of the multipole space will enhance 
the signal relative to the cosmic error. We use simple filter functions in / space to isolate 
different ranges of angular scales; low pass, Gaussian filters 



BiPS for different strength factors are shown in Figure HI As it has 
|30(], we can use different window functions in harmonic space, Wi, in 



W, 



G 



tG 



N"^ exp <^ - 



2/-I-1 
2/, + 1 



that cut power on scales {I > Ig) and band pass filters of the form 



wr 



2N- 



Jn 



2/ + 1 



exp 




(29) 



(30) 



2/, + 1 

where Jq is the ordinary Bessel function and A^*^ and A^'^ are normalization constants chosen 
such that, 



E 



(2/ + l)Wi 



21(1 + 1) 
i.e., unit rms for unit fiat band angular power spectrum Ci — mzri)- 
We use a simple x^ statistics to compare our BiPS results with zero. 



(31) 



2lT 



The probability of detecting a map with a given BiPS is then given by the probability 
distribution of the above x^- Probability of x^ versus the shear factor is plotted in Figure [5l 
Being on the conservative side we can constrain the shear to be (-^)o < 2.77 x 10^^° at a 
99% confidence level. 



12 



'n„-0'.5 arid (_£L)n 









", 


-0.5 and 


'I 


)^ -5.66e'- 


10 


4>;10 6 
















2x10 s 












































. 


T 






} 


I 


1 








^ ^ 


i 


















-2x108 

















'n„=0:5 and (ff^)„ =6,66e'-10 



-L--1- 













'n^ 


= 6,5 and 


(i.„^ 


l,le-10 




4xlO» 




















axio « 


- 
















- 





----- 


•---t- 


4- 


-:-- 




4-i- 




___L_L__ 


\-i' 


-2x10' 


- 



















FIG. 4: Bipolar power spectrum of Bianchi added CMB maps with strength factors of (-§')o = 
4.76 X 10"^°, 3.66 x 10"^°, 2.56 x 10"i°,1.25 x 10"^°. We have used W;^(4,13) to filter the maps 
in order to focus on the low-? (large angular scale) properties of the maps. 

V. CONCLUSION 

Bianchi templates, the characteristic temperature patterns in CMB, are imprinted by the 
unperturbed anisotropic expansion and have preferred directions which violate statistical 
isotropy of the CMB anisotropy maps. In this paper we study the consistency of the existence 
of a hidden Bianchi template in the WMAP data that has proposed in recent literature with 
the previous detection of null bipolar power spectrum in the WMAP first year maps |26| . 
We compute the bipolar power spectrum for low density Bianchi VII/j models embedded in 
background CMB anisotropy maps with the power spectrum that best fits the first year data 
of WMAP. We find non-zero bipolar power spectrum for models with (-§-)o < 2.77 x 10^^°. 
This is inconsistent with the null bipolar power spectrum of the WMAP data at (99%CL). 
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FIG. 5: Probability distribution of x^ versus the shear for Qq = 0.5 Bianchi template. Results 
of two different window functions, Wi^{25) (blue line) and VFj'^(4, 13) (red, solid line) are shown. 
By testing statistical isotropy of Bianchi embedded CMB maps, we can put a limit of {■^)o < 
2.77 X 10"^°(99%CL) on Bianchi VII/^ models. 

We conclude that correcting the WMAP first year maps for the Bianchi template may 
make some anomalies in the WMAP data vanish but this will be done at the expense of 
introducing other anomalies such as preferred directions and violation of statistical isotropy 
into the Bianchi corrected maps. 
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[64] All Killing vector fields form a Lie algebra of the symmetric operations on the space, with 
structure constants C^^. See [59] for mathematical preliminaries and details. 



[65] We split the matrix qab into two parts: its volume and the distortion part, qab = e^" (e'^^' 



AB ■ 



Where the scalar a represents the volumetric expansion whilst (3 is symmetric, trace free 3x3 
matrix and hence e^^ is given by the series 






r=0 



[66] Bianchi-I patterns with dark energy have been studied by |2z]. Also [36] have considered 

Bianchi VII/j models with dark energy. 
[67] Simulation of background maps and the spherical harmonic expansion of the resultant map is 



done using HEALPix package 
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